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Abstract . 


An  extensive  series  of  calculations  are  made  to  compare 
an  improved  closure  of  the  Born-Green-Yvon  (BGY)  equation  for  the 
electric  double  layer  of  primitive  electrolytes  to  existing  Monte 
Carlo  simulation  and  other  theories,  such  as  the  Modified  Poisson- 
Boltzmann  (version  5)  (MPB5)  and  the  Hypernetted  Chain/Mean  Sphe¬ 
rical  Approximation  (HNC/MSA) ,and  its  recent  improvements. 

In  contrast  to  these  theories  the  BGY  equation  satisfies 
the  contact  theorem  always.  Furthermore  the  bulk  pair  correlation 
functions  used  are  the  most  accurate  available,  namely  the  HNC 
bulk  results. 

The  results  show  very  good  general  agreement  to  the  compu¬ 
ter  work  for  high  densities  and  surface  charges.  In  particular 
we  are  able  to  go  to  higher  densities  than  those  hitherto  publi¬ 
shed  in  the  literature.  For  low  density  2-2  and  2-1  electrolytes 
the  agreement -is  not  as  good  as  for  the  other  theories. 


'  .  *  >  •  •  ■  •  *  .  ’ 
»  *  V  .  *  .  '  v'  . 


,  ,  *  ^  , 
V'  v*.*  *  * 


Introduction ■ 


The  primitive  model  of  the  electrode  interface  consists  of 
a  perfectly  smooth  surface,  which  represents  the  metal,  and  a  col¬ 
lection  of  charged  hard  spheres  representing  the  ions. All  of  it  is 
immersed  in  a  dielectric  continuum. 

Clearly  this  model  is  not  an  accurate  representation  of  the 
real  electrode-  electrolyte  interface,  but  it  can  be  simulated  in 
computers,  and  therefore  it  is  possible  to  assess  the  accuracy  of 
approximation  schemes  for  it. 

( 1  ) 

The  early  experimental  work  has  used  the  theories  of  Gouy 
Chapman^  and  Stern^  with  reasonable  success,  in  spite  of  the 
fact  that  these  theories  neglect  completely  the  excluded  size  effects 

and  correlation  effects  between  the  ions.  Recently,  various  theories 

(  4-'  0  ) 

have  been  proposed  to  include  these  correlations  in  the  treatment. 

One  of  the  primary  goals  of  this  exercise  is  to  develop  an  e- 
quation  caflble  of  handling  very  dense  and  higly  coupled  systems,  such 
as  the  molten  salts  ancT  the  non-primitive  electrolyte  (mixture  of 
hard  ions  and  dipoles)  against  a  metal  electrode. 

In  this  regime  it  is  most  important  that  the  contact  relati- 

(11,12) 

on 

kBT  f'  Pi(0)  =  Pb  *  *Eo  '  3jf  (1) 

be  exactly  satisfied.  In  (1)  k  T  is  the  Boltzmann  thermal  factor, 

O 

P  is  the  bulk  pressure,  e  is  the  dielectric  constant  (  e  z  30  for 

D 

water  in  the  primitive  model,  and  8  =1  for  the  non-primitive  case), 
<^(0)  is  the  contact  density  of  species  i  (  which  could  be  either 

3. 

an  ion  or  solvent  molecule  in  the  non-primitive  case)  and  the  sum¬ 
mation  is  over  all  the  species;  E  is  the  applied  electric  field. 


All  species  are  assumed  to  be  hard  spheres. 

The  BGY  equation  and  any  of  its  closures  will  always  sati¬ 
sfy  the  contact  relation.  Foe  this  reason  we  have  chosen  this  e- 
quation  in  this  work^’1^.  In  this  paper  we  study  and  extend  an 

f 

improved  ansatz  for  the  closure  of  the  3GY  which  we  proposed  earlier' 

Here  the  pair  correlation  function  is  forced  to  satisfy  the  local 

(91415) 

electroneutrality  condition  ’  ’  and  to  be  positive  definite. 

We  compare  the  results  of  this  theory  to  the  extensive  work  of  Yal- 


(1< 


(17) 


leau  and  coworkers'  -IAd  Snook  et  al.""  .  We  find  that  the  agreement 
is  excellent  for  the  values  of  the  contact  potentials  for  the  more 
concentrated  solutions  .  The  density  and  potential  profiles  are  in 
general  good  agreement  also. 

With  our  method  we  have  been  able  to  get  results  for  very  con¬ 
centrated  solutions  (up  to  4  M  of  a  1-1  electrolyte) .  In  these  ca¬ 
ses  we  observe  a  marked  layering  of  solute  against  the  electrode 
wall.  For  the  dilute  ionic  solutions  (  0.05  M  for  the  2-2  and  2-1 
electrolytes,  0.1  M  for  the  1-1  electrolytes)  the  results  are  in 
general  inferior  to  the  MPB5 . 

In  Section  2  we  describe  our.  basic  equations.  In  Section  3 
we  give  the  results  for  the  1-1,  2-2  and  2-1, 1-2  electrolytes. 

A  brief  discussion  of  the  general  conclusions  is  given  at  the 

end . 

Sect.  2:  Basic  equations. 

We  consider  a  model  in  which  the  ions  are  hard  spheres  of 

diameter  a  and  number  density  o.  .  The  solvent  is  a  continuum 
i  i 

of  dielectric  constant  e  ,  and  the  electrode  is  flat,  smooth  and 


-1 


has  the  same  dielectric  constant  e  . 


For  this  model  we  define  the  one  and  the  two  particle  di¬ 
stribution  functions 


0  (?  >  •  ^(D 


•lJCTl  •r2)  ■  SijU-2> 


(2)  , 
(3)  , 


where  r  s  (x^.y^.z^)  is  the  position  of  particle  1,  z^  being 
its  distance  from  the  wall 


The  pair  correlation  function  is  defined  by 


©1  Cl  ,2)  =  ^(1)  Q  (2)  g  Cl. 2) 


(4) 


and  we  also  need* 

h  ( 1 , 2)  s  g  (1 .2)  -  1  (5) . 

These  functions  must  satisfy  the  asymptotic  conditions 


lim: 

Z  — >oO 


lcu  “  »1- 

(6)  . 

bulk,  . 

ij(l’2  =  gij  ri2 

(7)  , 

, ,  .bulk  ,  . 

1J(1’2)  ■  hlJ  (rl2> 

(8)  , 

where  r12s  |ri-”  r2  I  This  means  that  all  the  correlations  far 
away  from  the  electrode  mu3t  be  equal  to  those  of  the  homogeneous 
bulk . 


Our  basic  approach 
as  a  modification  of  the 
Consider  the  first 


is  to  write  the  3orn-Green-Yvon  equation 

( 10) 

Poisson  -Boltzmann  equation 
BGY ^ ® ^  equation 


“kgT  ^lng  tCl) 


2 

j 


•il2)gu(l-2) 


M  ( r 

1  ij  12 


for  z  >  o  M  /  2 


where  the  pair  potential  U^(r  )  has  two  contributions 

U  .  .  (  r  )  =  U°.(r  )  +  e  e  /tv 

J-J  12  ij  12  i  y  *  12 


(in'  , 


where  is  the  short  ranged  part  cf  cue  ion-ion  interaction 

Ve  introduce  the  local  g lec - rostatic  potential 


0(1) 


-E  z  *1  dr  2  q  (2)  /  e  r 

°  i  J  2  i  9y  12 


(id. 


Adding  and  subtracting  the  gradient  of  the  second  term  of  (11) 
into  (9),  we  ejfeilnate  the  divergence  of  the  right  hand  side  of  (11) 
Using  also  the  relation 


h  U?j(ri2)  ‘  -V  4  (rl2'  (Z12/ri2>  U2! ' 


where 


•u  *  ( *i  *  ®j)/2 


(  13) 


Z12  =  Zl"  Z2 


(U), 


and  d(x)  is  the  Dirac  delta  function,  we  get  after  some  simple 
transformations : 

-V  ri  ln 

’  e!  vx*U)  -  V  J 


tfe  assume  first  that  the  inhomogeneous  pair  correlation  function 
is  known.  Then  equations  (11)  and  (15)  rorm  a  closed  pair  of  e- 
quations  for  the  unknown  0(1)  and  8^(1)  A  computationally  con¬ 
venient  form  of  these  is  obtained  by  noting  that  (11)  is  really  the 


^d?2etC2)gtJU,2)  ilrl2-  »u )(*ia/rl2) 


The  functions  f^(z)  are  found  by  requiring  that  the  local 
e  Lee  troneutral  i  ty  is  satisfied^  ’ 


-  e 


I  e 

J 


.  1  d"r  q  . 

Jj  2 


.(2)  h.  .(1  ,2) 
J-J 


The  effect  of  the  first  term  in  eq.  (20)  is  to  eliminate  the  reg: 

of  negative  g  .(1,2)  which  otherwise  could  occur  in  the  r.eighbcu: 

^  (10) 
hcod  of  the  electrode  surface 


Sect.  Ill:  Results. 


The  calculations  were  performed  by  the  technique  described 
( 10) 

elsewhere  ,  integrating  the  Poisson  differential  equation  using 
a  more  refined  version  of  the  predictor-corrector  integration  rou¬ 
tine  than  in  ref.  (10),  which  enhances  significantly  the  convergen¬ 
ce  rate 

Typically  3  to  S'  iterations  are  needed  for  the  lower  den¬ 
sities.  For  higher  molarities  we  need  at  most  10  iterations  to  reach 
convergence . 

In  our  calculations  we  have  used  the  same  parameters  as  tho¬ 


se  in  the  Monte  Carlo  simulations  of  Torrie  and  Valleau 


(  16) 


a  =  4.25  A,  «=  78.4,  T=298  K  The  adimensional  surface  cnarge 

*  2 

is  o  =  q  a  /e  where  q  is  the  surface  charge  density  (usually 

®  2  s  2 
in  units  MC/  cm  )  The  plasma  parameter  f  =  e  /e  a  -  1.6809  and 

|3  -  1/k  T.  Similarly,  the  adimensional  potential  drop  is 

O 


V  =  0e# 

where  $  is  the  potential  drop  across  the  interface  (usually  in 
millivolts).  We  have  explicitily,  for  the  above  parameters, 


=  88.7  a  *  (  C/cm  ) 


*  =  25.7  V*  (mV) 


a)  1-1  electrolytes . 


We  have  computed  the  potential  and  density  profiles  for  so- 

* 

Lutions  ranging  from  0.1M  to  4M,  and  surface  charge  a  ranging 
from  0.1  to  0.7.  The  diffuse  layer  potential  drop  for  all  the  re¬ 
ported  cases  is  given  in  Table  I  where  we  also  display  the  Monte 
Carlo  results  and  those  of  competing  theories.  The  most  exter.si- 
ve  results  are  those  of  the  MPB5^  and  the  HNC , MSA^ 

It  is  noteworthy  to  say  chat  the  MPB5  uses  a  Debye-Huckel 
like  bulk  correlation  functions,  while  the  good  agreement  shown  by 

the  HNC/ MSA  is  obtained  through  the  use  of  bulk  correlation  func- 
■ 

tions  which  are.  poorer  than  those  used  in  the  HNC/ HNC  which  is  not 
a  very  good  theory  (  ref.  6,  1979). 

(19) 

In  our  work  we  use  most  accurate  bulk  correlation  functions 
which  are  in  good  agreement  with  earlier  reported  work  of  Rasaiah 
and  Friedman^20  ,  for  those  cases  that  these  last  authors  have  stu¬ 
died  . 

Furthermore,  our  work  will  always  satisfy  the  contact  rela¬ 
tion  (1)  within  the  numerical  accuracy  of  the  calculation .which  is 
generally  less  than  196  in  all  cases  investigated  and  typically  clo¬ 
se  to  0.396.  This  is  particularly  important  for  high  densities  and 
Low  couplings  .  Again,  neither  the  MFB5  nor  the  HNC/ MSA  will  sati¬ 
sfy  this  relation  exactly.  For  the  0.1M  the  results  are  given  for 

*  *  as 

a  0  1,0  2  and  0.3  :  the  potential  drop  V  are  not  accurate  as 

those  of  the  MPB5 . 

In  figs. 1-2  we  show  the  potential  and  density  profiles  for 

* 

the  case  0  =0.3.  The  agreement  of  the  BGY  is  better  than  that  of 

the  MGC ,  and  in  general  quite  good.  We  do  not  have  data  for  this 
case  from  the  MPB5.  In  the  low  density  regime  we  experience  nume¬ 
rical  difficulties  with  the  electroneutrality  condition  (22),  which 
could  be  of  technical  nature. 


se  that  the  overall  agreement  with  the  Monte  Carlo  is  quite  good. 

The  BGY  is  slightly  better  than  the  HNC ,  MSA  and  slightly  worse  tna 

*  * 

the  MF35  for  surface  charges  below  0  =0.5. For  O  =0.55,0.6 

the  agreement  to  i»IC  is  good  (  see  fig. 3) 

( 21 ) 

Recently  Nielaba  and  Forstmann  have  reported  a  new  appro 

ximation,  the  HNC/LMSA,  where  the  correlation  functions  are  compu¬ 
ted  using  the  MSA  at  a  local  average  concentration.  The  results  ;• 
the  potential  as  a  function  of  the  reduced  charge  are  s  icwn  in  fig 
3,  where  it  is  also  visibLe  the  failure  of  the  HNC, MSA  for  high 
charge  densities. 

» 

In  fig.  4  we  show  a  comparison  of  the  case  with  a  =0.141 
(1-1  case  at  1.0M)  with  MC  results.  The  agreement  is  very  good  and 
of  comparable  accuracy  than  MPB4^16^. 

For  highest  surface  charge  we  show  both  the  potential  profi¬ 
le  (  fig.  5)  and  the  charge  profilesC fig.6) .  It  is  interesting  to 
note  that  the  BGY  reproduces  quite  well  the  qualitative  features 
of  the  counterion  density  profile,  in  particular  the  dens  ity  oscil 

lations  near  z-l .  The  agreement  is  however  not  as  good  as  that  of 

( 21 ) 

Nielaba  and  Forstmann  ,  but  we  do  not  have  any  adjustable  para¬ 
meter  in  our  calculation. 

This  feature  of  density  oscillation  due  to  a  layering  near 
the  electrode -surface  is  also  present  in  the  higher  density  regime 
For  the  2.0M  case,  already  investigated  in  ref  .(13),  the  agreemen 
in  the  potential  drop  is  quantitative,  while  the  density  profiles 
are  in  remarkably  good  agreement  with  the  MC  results. 

In  fig.  7  we  show  the  f.  ntial  drop  as  a  function  of  the  su 
face  charge  for  2.0M  and  3 .0  M  .  On  the  basis  of  these  results  we 
also  are  able  to  compute  the  differential  capacitance  of  the  dif- 


fuse  layer,  a  quantity  of  interest  to  the  electrochemists. 

In  fig.  8  we  show  the  potential  and  density  profiles  for 
3.0  and  4.0  M  solutions  These  show  very  pronounced  layering. 

The  potential  profile  also  shows  a  minimum  which  is  deeper  for 
the  higher  density  case. 

b )  2-2  electrolytes. 

We  have  computed  results  for  two  concentrations:  0.05  M 
and  0.5  M.  The  results  for  the  potential  drop  are  displayed  in 
Table  II.  Here  again,  the  results  for  the  more  dilute  case  are  not 
as  good  as  those  of  the  MPB5 ,  while  the  0.5M  case  is  in  excellent 

i 

agreement.  We  experience  also  here  severe  problems  with  the  elec- 

a 

troneutrality  relation,  which  could  again  be  of  technical  nature. 

The  results  for  the  potential  drop  are  shown  in  fig.  9,  whe¬ 
re  we  also  display  the  differential  capacitance. 

In  fig.  10  we  show  the  charge  distribution  of  the  0.05M  (2-2'! 
case.  In  spite  of  the  poor  agreement  of  the  potential  drop,  the 
density  profile  are  in  qualitative  agreement  with  the  MC  results. 
The  counterion  profile  is  displaced  to  the  right,  and  the  coion 
profile  has  a  peak  at  the  correct  place. 

For  the~0.5M  case  the  results  for  the  potential  and  density 

profiles  are  shown  in  figs.  11  and  12  The  agreement  here  is  good 

In  particular, 

both  in  the  potential  drop  and  in  the  density  profiles,  the  poten¬ 
tial  profile  is  in  good  agreement  near  the  wall,  but  the  3GY  is 
higher  than  the  MC  far  from  the  wall  The  MPB5  is  always  lower  than 
the  MC  result. 

c )  Asymmetric  case. 


Two  concentrations  were  investigated :0 .05M  and  0.5  M  for  the 


2-1  and  1-2  cases.  The  results  tor  tne  potential  drop  are  displayed 
in  table  III.  Again, here  the  agreement  for  the  more  dilute  case  is 
poorer,  presumably  due  to  technical  reasos, 

Fig.  13  shows  the  results  of  the  potential  drop , as  well  as 
chose  of  the  differential  capac i tance , for  the  asymmetric  case.  As 
can  be  seen  the  agreement  is  better  for  singly  charged  counterions. 

Figs.  14  and  15  show  the  potential  and  density  profiles  for  the  Q.5M 

♦ 

(1-2,  i.e.  doubly  charged  counterion)  at  surface  charge  a  =0.1704. 
While  the  potential  profile  is  very  similar  to  that  of  the  MP35, 
the  density  profiles  are  slightly  better  and  in  excellent  agreement 
with  the  MC  profile. 

« 

We  have  also  performed  the  calculation  of  the  0  -0 ,  0.5  M 
case,  recently  reported  by  Torrie,  Valleau  and  Outhwaite^ .  Fig. 

16  shows  the  potential  profile,  which  is  slightly  inferior  to  the 
MPB5  result,  and  fig.  17  the  density  profiles  ,  which  are  slightly 
better  than  the  '!P95 


Concluding  remarks . 


In  the  present  work  we  have  presented  an  extensive  discussion 

of  the  use  of  the  BGY  for  charged  interfaces.  The  BGY  equation  is 

an  exact  relation  but  requires  a  closure,  which  is  an  ansatz  for  the 

inhomogeneous  pair  distribution  function  g.  .(1,2). 

1 J ( i 3) 

In  this  work  we  have  used  a  construct  that  removes  some 

of  the  unphysical  regions  with  negative  g  (1,2)  that  occurred  when 

(9,10)  ^ 

the  modified  superposition  ansatz  ’  was  used.  The  results  are 
in  very  good  agreement  for  high  salt  concentrations  but  definitely 
poorer  for  dilute  solutions.  The  addition  of  a  repulsive  part  to 
the  pair  interaction  near  the  wall  takes  into  account  higher  repul¬ 


sive  interactions  due  to  a  crowding  effect.  We  must  remark  that  the 


affect  of  this  term  is  in  general  not  very  large.  However,  it  in¬ 
creases  very  significantly  the  stability  of  the  iteration  prccedu 
re  . 

The  poor  performance  for  dilute  solutions  is  attributed  to 
difficulties  with  the  electroneutrality  rule,  which  could  also  be 
of  a  technical  nature.  We  will  investigate  this  point  in  the  futu 

The  main  objective  of  our  work  is  to  find  an  equation  that 
will  perform  well  in  the  high  density-high  coupling  regime  whi  ;h 
is  that  of  the  molecular  solutions  or  molten  salts  near  a  metal 
electrode.  Our  equation  is  indeed  well  behaved  in  this  regime, and 
we  have  been  able  to  obtain  physically  reasonable  solutions  at  3. 
and  4.0M,for  fairly  high  surface  charges. 
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Electrolyte  0.05M 


0  170-1  0.45  0.460  0.487  0.574  5  10.65  10.787 


Table  I:  (a)  Diffuse  layer  potential  drop  yp  (0)  for  1-1  elec¬ 
trolytes.  MC  results  ;G.M.Torrie  and  J.P.Valleau  (1980) 
(ref.  16) ;MPB5 :C .W.Outhwaite  and  L.3.3huiyan  (1983), (re 
8)  :HNC/MSA:M.  Lozada-Cassou  e_t  al  .  (  1982), (ref.  7) 

(b)  As  calculated  from  the  r.h.s.  of  eq.  (1). 

♦ 

Table  II:Diffuse  layer  potential  drop  yf  (0)  for  2-2  electrolyte 
MC  results :G  M.Torrie  and  J.P.Valleau  (1982),  (ref.  16): 
MPB5  and  HNC/MSA  same  captions  as  Table  I.  The  HNC/MSA 
results  at  005  M  have  been  taken  from  a  graph. 

* 

Table  III:  Diffuse  layer  potential  drop  yp  (0)  for  1-2  and  2-1 
electrolytes.  Same  captions  as  Table  II. 


an  i 


"ig.  : 

rig.  2 

Fig.  3 


Fig.  4 

Fig .  5 
Fig  6 

Fig.  7 

Fig.  8 


Fig.  9 


Potential  profile  for  a  1-1  electrolyte  at  cO-lM 

* 

a  =0.3.  Dots,  MC  results(ref . 16) :  full  Line,  this  work. 

Density  profiles  for  the  case  of  fig.  1.  Dots  and  crosse 
MC  results; full  line,  this  work. 

* 

Diffuse  layer  potential  drop  y>  (0)  (  0/2  is  taken  as  cc 

dinate  origin)  as  a  function  of  the  surface  charge  der.si 

* 

ty  a  .  Dots,  MC;dashed  line,  HNC,  LMSA;dashed-doubie  dot 
1  ine  ,MPB5  ;  dash-dot  1  ir.e  ,  HNC/ MSA ;  full  line,  this  work. 
Also  shown  is  the  differential  capacitance  in  the  3GY(oh 
work),  whose  scale  is  on  the  rightCin  8.854  fiF /  crr.~  units 

Density  profiles  for  the  1-1  electrolyte  at  csl.J.-i  and  a 
=0.141.  Dots  and  cresses, MC  results;full  line,  dashed  li 
ne ,  this  work. 

# 

Same  captions  as  fig.  1  except  that  C=1 .OM  and  a  =0.7. 

Density  profiles  for  fig.  5  case.  Dots  and  triangles,  MC 
resultsjopen  circles  and  crosses , HNC/ LMSA;  full  line,  da 

shed  line,  this  work. 

*  * 

V  (0)  vs.  a  for  the  1-1.  case  at  c=2.0M  and  3.0M.3GY  re 
sults:fulx  line,  2.0M;dashed  line,  3.0  M.  Also  shewn 
are  the  differential  capacitances( same  units  as  fig  3). 

Potential  and  density  profiles  for  3.0  and  4.0M  con¬ 
centration  of  the  1-1  electrolyte.  Frop  top  to  bottom: 

potential,  counterion,  coion  profile.  Dash-dot  line,  3. OH 

*  * 

at  O  =0.485. Full  line,3.0M  at  a  =0 . 396 ; dashed  Line,  4.0 
M  at  a* =0.56. 

•  * 

ip  (0)  and  diff.  capac.  vs.  a  for  the  2-2  electrolyte 

at  c=0.5M.(in  the  same  units  of  fig. 3) 


0 


=0.284.  Dots  ar.d  tresses,  M  ’  reui  :s .  F  ail  Hne.  ...  r 

dashed  line,  Modified  Gcuy- Chapman . 

Fig.  11  :  Potential  profile  for  the  2-2  case  at  c=0.5M  and  a  =0.1' 

Dots,  MC  results;  full  line,  this  work  ;  dashed  line,  .IF 

Fig.  12  :  Density  profiles  for  the  case  of  fig.  11.  Dots  and  cres¬ 

ses,  MC  results;  the  rest,  same  captions. 


rig. 
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CO)  vs.  a  for  1-2, 2-1  electrolytes  at  different  mo  la 

* 

rities.  Counterion  single  charged  for  a  >  0. doubly  .narg 

* 

for  tf  <  0  .Concentration  c=0.5M:dots,  MC  results  :  full  Hr. 
this  work.  c=0.05M:  crosses,  MC  results ; dashed  line,  thi 
work.  Dash-dot  line:  differential  capacitance  for  the 
c=0.05M  case  (same  units  of  fig  3). 


Fig.  14 

I 


:  Potential  profile  for  the  l-2(doubly  charged  counterion) 

* 

case  at  c=0.5M  and  a  =0.1704.  Dots.MC  results;  dashed  . 1 
MP85;  full  line,  this  work. 


Fig.  15 


Fig.  16 


Density  profiles  for  the  case  of  fig.  14.  Dots  and  cross 
MC  results.  For  the  rest  same  captions. 

* 

Potential  profile  for  the  1-2  case  at  a  =0  and  c=0.5M. 
Dots,  MC  results;  crosses,  MPB5 ;  full  line,  this  work. 


Fig.  17  :  Density  profiles  for  the  case  of  f ig . 16 . S ingle  charged 

ionf  crosses,  MC  results ; dash-double  dot  line,  MF35; 
dashed  line,  this  work. Doubly  charged  ion:  dots,  MC  resu 
dash-dot  line,  MPB5;full  line,  this  work. 
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